Magnetism and Charge ordering in TMTTF2-PF6 organic crystals 
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Using a combination of Density Functional Theory, mean-field analysis and exact diagonaliza- 
tion calculations we reveal the emergence of a dimerized charge ordered state in TMTTF2-PF6 
organic crystal. The interplay between charge and spin order leads to a rich phase diagram. Coex- 
istence of charge ordering with a structural dimerization results in a ferroelectric phase, which has 
been observed experimentally. The tendency to the dimerization is magnetically driven revealing 
TMTTF2-PF6 as a multiferroic material. 
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In recent years, organic molecular crystals have 
emerged as a fascinating class of materials with immense 
potential for applications. Structures in these materials 
are held together by van der Waals and hydrogen bonding 
as opposed to much stronger covalent and ionic bonding 
in conventional crystals. As a consequence, these systems 
present complex phases and phenomena which still re- 
quire a microscopic understanding, ranging from anoma- 
lous positive and negative thermal expansion[l| to high- 
temperature ferroelectricity with large polarization[2] 
and even superconductivity 0. Moreover, their light, 
flexible and nontoxic character makes them ideal can- 
didates for future applications. 

The TMTTF2-X series of charge transfer salts shows 
a very rich phase diagram when either hydrostatic or 
chemical pressure is applied 0, [Hj with electronic phases 
ranging from spin/charge density waves to a spin-Peierls 
state. The case of TMTTFa-PFg is pecuhar. At ambi- 
ent pressure it is a Mott insulator that develops a charge 
ordered (CO) state with a 4ki? charge localization and 
the CO transition coincides with the onset of ferroelec- 
tric order Moreover by lowering the temperature 
this state can be turned into an antiferromagnetic (AFM) 
state or a Spin-Peierls phase. 

In this paper we use a variety of theoretical and com- 
putational methods - namely, density functional the- 
ory (DFT), mean-field (MF) and exact diagonalization 
(ED) - to describe the low-temperature phase diagram 
and the emergence of a ferroelectric state in TMTTF2- 
PFg. We show that the interplay between charge- 
ordering, spin-ordering and structural dimerization leads 
to a rich phase diagram including a dimerized charge- 
ordered phase which breaks inversion symmetry and con- 
sequently is ferroelectric. 

We begin by presenting the results of DFT calcula- 
tions, which are carried out using the generalized gra- 
dient approximation to the exchange-correlation accord- 



ing to Perdew-Becke-Ernzerhof (PBE) [8|. The electronic 
structure of TMTTF2-PF6 is computed using Quantum 
Espresso package Q and a tight-binding representation of 
the bandstructure is built using WannierQO to com- 
pute the maximally localized Wannier orbitals. The cal- 
culations for the low ternperature centrosymmetric struc- 
ture of TMTTF2-PF6 [uj provide a simple picture: the 
HOMOs (highest occupied molecular orbitals) of the two 
TMTTF molecules hybridize, giving rise to bonding and 
antibonding bands around the Fermi level, with charge 
transfer to the PFg anions. The ratio 2:1 between cations 
(TMTTF molecules) and anions (PFg group) in the unit 
cell leads to a commensurate band filling of 3/4 of the 
TMTTF band. The system is metallic within PBE with 
a bandwidth W of ^ 1.0 eV and it shows a quasi one- 
dimensional character along the stacking chain direction 
a as reflected in the band structure. The inability of 
PBE to reproduce the experimental insulating behavior 
is a very well known shortcoming of this approach when 
strong local and short-ranged correlations are relevant. 

In order to include such correlation effects we model 
TMTTF2-PF6 by a two-dimensional extended Hubbard 
model with the Hamiltonian: 
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Here, c^^ and cj^. are annihilation and creation opera- 
tors for electrons at site i with spin a, tij denote the bare 
hopping strengths and Ui are the distortions along the 
stacking chain direction a of the TMTTF ions from their 
equilibrium position. The electron-lattice coupling pa- 
rameter is a, which favors dimerization. U is the strength 
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FIG. 1. (Color online) a) PBE band structure of TMTTF2- 
PFe in the low temperature centrosymmetric phase. The 
high symmetry directions are labeled as: F = (0,0, 0), X = 
(i,0,0), Af = (i, |,0), 7? = (0, i,0). b) Schematic structure 
of various hopping parameters in TMTTF2-PF6 molecular 
crystals in the a,b plane along a {ta,t a) and b (ti,), and the 
diagonal directions {ta+b,t a+b)- 



^) V'(eV) U=2.2(eV) 



0.4 - 



b) 



0.2 - 



0.0 - 




I I I I I I I I I I I 
0.0 0.2 0.4 0.6 0.8 1.0 

V(eV) 



Jffft 

fttt 
, fttt 



FIG. 2. (Color online) (a) Phase diagram showing the mag- 
netic and charge ordered phases in the parameter space of V 
and V' (eV). M and I label Metallic and Insulating ground 
states respectively. A schematic view of charge and magnetic 
orders appearing in the phase diagram is shown in (b) for a 
unit cell 4x2. 



of the on-site Hubbard interaction, V is the inter-site re- 
pulsion strength between TMTTF sites along a, and V 
is the inter-site repulsion strengths along b and along the 
diagonal a + b (see Fig. [TJ))). 

The hopping parameters are extracted mapping the 
corresponding tight-binding model on the basis of max- 
imally localized Wannier orbitals from PBE band struc- 
ture. The only non- negligible hoppings are shown in 
Fig. \Sp) and their values are ta^O.22 eV, t'a=0.20 eV, 
tb=OM eV, ta+b=-OM eV, t a+b^-Om eV, confirming 
the high anisotropy of the electronic properties. These 



values agree well with experimental estimates of Ref. [12 
and are rather larger than those estimated on the ba- 
sis of Hiickel calculations [ll| in which the effect of the 
anions and of the screening of the crystal are not taken 
into account. The slight difference between the two hop- 
ping parameters along the a axis ta and t a is due to the 
crystal arrangement of PFg anions which determines a 
dimerization of the TMTTF molecules (the two inequiv- 
alent distances S and L in Fig. [T] have difference of ~ 
0.04 A[ii[). In terms of our model ([T]), this corresponds 
to a ^1.2 assuming ta—t a=0.21 eV. 

The Hubbard U is the screened energy cost required 
to doubly occupy a molecular site with holes. This en- 
ergy cost can be evaluated explicitly from the total en- 
ergy of an ion with different number of electrons and 
the polarization energy [l^l- We calculate the variation 
of the eigenvalues and their polarization energy of the 
(hole) charged molecules placed inside a cavity of an ho- 
mogeneous dielectric medium with dielectric constant of 



PBE functional. 

Since little is known theoretically about the inter-site 
interaction parameters V,V' we adopt a phenomenolog- 
ical approach studying the phase diagram letting these 
parameters vary in the typical range of other molecu- 
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using the GAMESS package 



1. We estimate ;7=2.2 eV 
15| using TZV basis set and 



lar crystals [16[. An experimental estimate V c^i 2t [17 1 
would imply V ~ 0.4 eV, which lies in the middle of the 
range we consider. 

Employing a mean-field Hartree-Fock (MF) decompo- 
sition for the Hubbard U and V,V' terms, we solve the 
model in Eq. [T]with the hopping parameters listed above 
self-consistently for different initial conditions and retain 
the lowest-energy solutions using a unit cell of 4x4 sites. 
In Fig. [21 we show the groundstate phase diagram in the 
parameter space of V and V . We obtain a variety of 
phases that emerge from the interplay between the on- 
site and inter-site repulsions. The main tendency at MF 
level is the charge ordering which indeed occurs at any 
finite values of V and V . We find essentially two struc- 
tures: (i) the CO pattern with wave vector q= (tt, 0) 
(phases 6, 7 in Fig. [2]), stabilized with increasing V at 
finite V' , which can be directly connected with the fer- 
roelectric state for X=PF6, and (ii) the CO pattern with 
wave vector q= (tt, tt) (phases 2, 3, 4, 5 in Fig. [5]), stabi- 
lized with increasing V' at V=0, which is similar to the 
antiferroelectric state observed with X=SCN [ISj . 

To a large extent, the spin-ordering follows the CO 
pattern and combines with it leading to the rich phase 
diagram we have drawn. For small values of V the fer- 
romagnetic (FM) order it-fr-lt-tT along the stacking chain 
direction a is the groundstate (phase 1 in Fig. [5]) with 
a lower energy with respect to the AFM state lt-'fr--lj---ll- 
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which is realized without CO at higher energy. The FM 
solution can be tuned into a CO checkerboard state at 
finite V and V'=0 (phase 2 in Fig. [5]). Increasing V at 
finite V an AFM state with charge and spin modulation 
along direction a establishes. In such state the spin pat- 
tern is 'fr-t-JJ--4 (where the single arrows refer to a small 
magnetic moment and the double arrows to a large one) 
and staggered charge pattern 6'^-S~-S~^-S~ {S^ indicat- 
ing a positive or negative charge imbalance with respect 
to the average) leading to phases 3 and 6. At larger V the 
AFM state turns into the extreme picture '^-O-ij.-O with 
alternating singly and doubly occupied TMTTF sites, the 
latter having zero magnetic momentum (phases 4, 5 and 
7) . It is important to realize that in our two-dimensional 
model the behavior along the a axis combines with the 
CO pattern along the b direction that we discussed above 
leading to five different phases (numbered from 3 to 7). 
All these phases are characterized by 2kp spin ordering 
and 4kp CO along the a-direction. They are consistent 
with the experimental evidence of spin and charge [1] 
ordering in TMTTF salts and well agree also with calcu- 
lations in the one-dimensional extended Hubbard model 
[20| . Our calculations confirm these results in a more 
realistic two-dimensional model, and show that it is of 
great importance to take into account the complicated 
details of the interchain couplings to study the compe- 
tition of the different phases. The degree of anisotropy 
of the intersite interactions in TMTTF2-PF6 crystals is 
more complicated than the one used in our model how- 
ever our approach captures both CO patterns revealed 
from experiments in this class of materials. 



Lattice effects can accompany the CO phase [21| and 
indeed a Akp bond dimerization has been reported ex- 
perimentally for the CO phase of TMTTFa-PFg [ll|. 



Atomic displacements, which break the inversion symme- 
try, can give rise to ferroelectricity once they are coupled 
with a coherent charge order as in the case of TMTTF2- 
FFg. Here in the following we discuss, as a natural gener- 
alization of the dimerized hopping structure due to the in- 
equivalent distribution of PFg anions assumed so far, the 
possibility of formation of dimers between TMTTF sites 
along the stacking chain direction a as consequence of 
minimization of the exchange energy: the onset of AFM 
ordering pins the dimer formation tendency. For the sake 
of definiteness, we assume in the following V' = (ob- 
viously V is expected to be smaller than V because of 
the larger distance involved). 

We now turn to the model of Eq. [T] considering a sym- 
metrized hopping (ta = t a=0.21 eV) and we vary a from 
to 2 in order to explore the interplay between CO and 
dimerization. The energy is minimized as a function of 
the classical distortions Ui. In Fig. [3] the phase diagram 
shows a transition from an undistorted state (phases 6' 
and 7') to a dimerized state (phases 6 and 7) with in- 
creasing a. The key point is that a larger a is needed to 
stabilize dimerization when CO is stronger. The value of 
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FIG. 3. (Color online) Phase diagram showing the dimeriza- 
tion of magnetic and charge order phases in the parameter 
space of a and V. Here we have used V' — V/2. 



a required to favor the transition to a dimerized state is 
in good agreement with our DFT estimate in the whole 
range of V^s. Our results complement previous similar 



reports 22|-|24| within a realistic framework in which the 
bandstructure is derived from ab-initio calculations. 

The 4ki?-CDW of a Wigner type with dimerization of 
TMTTF lacks inversion symmetry and is therefore fer- 
roelectric as experimentally reported Q . Equivalent dis- 
placements (As) of molecular sites in the formation of 
dimers along the stacking direction with disproportiona- 
tion of charge {S~^-S~) create local dipoles of value 25 As 
at each dimer, which sum up along the lattice in the a 
direction resulting in a finite polarization in the ferroelec- 
tric CO at (7r,0). In this scenario the anion PFq~ is ex- 
pected to move closer to the molecular site with less elec- 
trons ((5+) making the CO stronger [25|- The CO states 
(phases 6, 7 in Fig. [2|) are characterized by a unit cell con- 
taining an essentially singly occupied TMTTF molecule 
(A) and an essentially doubly occupied one (B). In these 
states the spin-minority states at site A are empty and 
can be occupied by minority electrons jumping from site 
B creating exchange striction of the same kind found in 
novel organic multiferroic materials [26| . The stronger 
the magnetic coupling, the easier it is for the system 
to dimerize. Therefore the magnetism is controlling the 
ferroelectricity. The symmetry of the insulating state 
is such that magnetic ordering pushes a charge-ordering 
pattern from site-centered towards bond-centered [27 1. 



In this intuitive picture of the origin of the ferroic state 
we treat electrons as localized however the covalency of 
the electronic state pumps the charge through the crys- 
tal by quantum effects |28| . It is interesting to observe 
that moderate values of '~ 0.5eV are sufficient to lead 
to this extreme charge disproportion despite the larger 
value of U. 

Here we remark the importance of using a full two- 
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FIG. 4. (Color online) (a) Spin and charge structure factors 
S{q) and C{q) at some relevant wave vectors (The two S{q) 
for different momenta lie on top of each other), (b) Charge 
gap Ac at different values of V and V' = V/2 in MF and ED 
schemes. 



dimensional modeling to properly account for the mag- 
netic and charge properties. A one-dimensional de- 
scription would indeed miss the competition of anti- 
ferroelectric phases, whose occurrence depends on com- 
plicated details of the interchain coupling. The present 
MF scheme however does not take into account quan- 
tum fluctuations, and hence overestimates the stability 
of magnetic and charge orderings. In order to take into 
account such fluctuations, we solve ([IJ by Lanczos ED 
on a 3/4- filled finite square cluster of 4x4 sites (N) with 
the same parameters of our MF analysis. Within this ap- 
proach we calculate the charge gap (Ac) and the struc- 
ture factors jfJ2i^^Pi'^1^-i)i''^if ^ "-ij.)) where the plus 
and minus sign correspond to the c harg e {C{q)) and spin 
{S{q)) structure factor respectively [29|. 

In Fig. m we show C{q) and S{q) for the momenta 
that maximize the response functions as a function of V. 
The ED C{q) confirms indeed that V favors CO with 
wave vector q=(7r,0) (^kp) for V > 0.2 eV (see Fig. gji)), 
while for small V it is competitive with the q=(7r,7r) CO 
(not shown). S{q) displays a smaller enhancement as a 
function of V and two magnetic orderings are empha- 
sized: (i) an AFM state with wave vector q=(^,7r) {f- 
it-JJ'-JJ-) [30j and (ii) another AFM state with wave vec- 
tor q=(|,|). The two magnetic orderings are consistent 
with the one we find in MF for X=PF6 and observed ex- 
perimentally for X^SbFg and X=SCN jJl] respectively. 
Coupling charge and magnetic states we can recover the 
same phases 6 and 7 obtained in MF. The charge gap 
Ac is calculated averaging over twisted boundary con- 
ditions [sij in order to reduce finite-size effects and the 
MF values are shown in Fig. IDd) for selected values of 
V. The charge gap Ac is finite over the entire range of 



V. Increasing the lattice size is expected to reduce the 
value of the charge gap moving toward the experimental 
value [1^ and suggesting that a small value of V may 
be relevant for such organic structures. Moreover the 
temperature effects not taken into account in the present 
study soften both spin and charge orders. 

The tendency toward formation of dimers along the 
stacking direction a is qualitatively maintained within 
ED calculations. The competition between two insulat- 
ing phases may lead to a metallic state in the V-V phase 
diagram which results from the geometrical frustration 
of the long-range Coulomb interactions [sll but due to 
the anisotropy and large on-site Coulomb repulsion the 
TMTTF2-PF6 salt is in a dimer Mott insulating state. 
In conclusions with a combination of different techniques 
we have highlighted a peculiar novel interplay between 
charge, spin and lattice degrees of freedom in the low tem- 
perature phase diagram of TMTTF2-PF6 crystals. The 
magnetic and charge orderings coming out from MF and 
ED analysis are consistent and agree with experiments 
suggesting that the competition between local and non- 
local Coulomb repulsion is the key to understand the 
properties of this compound. The charge ordered state 
in conjunction with the spin ordering makes the system a 
natural candidate for being a new multiferroic material. 
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